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ABSTRACT 


De^ntng  and  testing  dam^e-resistant  Navy  lequirea  meaanieniNit 
oS  abscdate  ship  motion  caused  by  attacks  from  ndties,  d^th  diaiges,  totpe^ies. 
etc.  This  rq>ort  ^es  a  mefliod  of  designing  flectronic  ist^iation  systems, 
stabilized  by  feedback,  which  may  be  used  to  obt^  abstrfote  idiip  motion  by 
integrating  accderometer  oatputs.  Design  data  are  given  for  sevod  variations  of 
die  genoal  integration  mediod  presented,  and  an  example  of  eadi  variation  is 
demonstrated  by  obtaining  di^lacemait  from  accderation  by  dooUe  int^ntioii. 

ADMINISTRATIVE  INFORMATION 

The  woHc  presented  in  this  tepori  was  began  under  Defense  Atinnic  Support 
Agency  Project  Number  K-llBAXN,  Task  Area  XSOl  and  was  completed  under  Navy 
Picjcct  Number  S-F3S.422.110,  Task  15050. 

INTRODUCTION 


BACKGROUND 

Navy  ships  and  equipment  must  be  designed  to  resist  damage  frmn  underwater 
explosions  due  to  mines,  depdi  charges,  torpedoes,  etc.  The  design  of  damage-resistant 
diips  and  equipment  which  meet  Navy  shock  and  vibration  specifications  requires  know¬ 
ledge  of  the  relatioQship  between  underwater  explodons  and  damage.  To  obtain  this 
knowledge,  dynamic  measurements  of  such  mechanical  parameters  as  dis{dacement,  vdocity 
and  accderatioi)  must  be  made  during  underwater  explosion  experiment  on  sh^  Hie 
severe  shock  mivironment  in  which  sensing  devices  must  operate  and  tfie  lack  of  con"-m- 
ient  stationary  refaences  for  measurements  of  absdute  motion  durii^  the  expoinient 
create  uniqae  instrumentation  proUems  in  underwater  explosions  work.  This  reptwt  deals 
with  the  prot^em  of  absolute  motion  measurement  during  explrsion  experiments  and 
presents  a  solution  to  critical  aspects  of  the  problem. 

During  past  underwater  explosion  tests,  vdocity  and  displacement  measurements 
at  various  location!)  on  the  target  ship  have  been  made  inertially  by  seismically  conectii^* 

I.  Wsftei,  lOR.,  **A  ProcedBR  for  Uic  Cotiectioa  of  Velocity  Meter  Rcccndt,"  Darid  Tmlor  Model  Baia  Reeott  i960 
iJiin  1965>. 

A  comfdele  listing  of  references  s  gtrea  on  page  SI . 
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and  integrath^  the  oatpat  of  vdocity  meters^  and  photographicaOy  by  usii^  fixed  refer* 
OKcs.  The  limited  di^dacement  rai^e  of  vdocity  meteis  and  the  i»oMeni  of  providiiig 
a  stationary  camoa  and  referotce  for  {diotographic  measurement  have  shown  a  need  for 
devdoping  better  methods  for  making  these  motion  measmemaits. 

Until  recently  accderometers  have  not  been  availaUe  that  have  diaractmistics 
saitaide  for  accuratdy  senang  diip  motions.  However  recent  improvements  in  avaflaUe 
accderometers  allow  their  use  in  many  situations.  These  accdertMneters  provide  a  v<dt^ 
ou^t  prr^rtional  to  die  sensed  accderation.  Vdocity  and  displacement  may  be  obtained 
fiom  the  measured  accderation  by  integration  and  doulde  integration,  reqiecdvdy.  For 
the  time  duration  required  in  underwater  exnlosions  work,  the  dynamic  ranges  of  most 
general  purpose  instrumentation  tape  recorders  are  not  large  enou^  to  record  accderation 
mth  the  accuracy  needed  for  accurate  integration  and  doulde  integration.  Therefore  the 
integration  {nrocess  used  must  operate  directly  on  the  signal  generated.  Snce  die  accder* 
ometer  lus  a  voltage  output,  it  may  be  dire<'tly  integrated  using  dectronk  integrators. 
Vdocity  and  diiplacement  initial  values  and  the  accderation  component  due  to  gravity 
are  unknown  because  target  motion  prior  to  die  test  is  present  Therefore  open  loop 
operation  of  the  integration  system,  which  requires  starting  the  integrators  with  correct 
initial  values  immediately  before  the  test  is  not  practical.  A  pracijcal  integration  system 
must  be  capable  of  being  started  sometime  ahead  of  the  test  so  as  to  automaticaUy 
determine  vdodty  and  displacement  initial  values  at  the  beginning  of  the  transient 
accderation  integrated. 

OBJECTIVE 

The  purpose  of  this  rqiort  is  to  present  and  demonstrate  a  practicable  method 
of  performing  multiple  integration  of  a  transducer  output  signal  that  will  be  applicalde 
under  field-test  conditions.  To  meet  this  objective,  design  data  are  given  for  several 
varktions  of  die  general  integration  method  presented,  and  an  example  of  each  variation 
is  demonstrated  by  obtaining  di^lacement  from  acceleration  by  dooUe  inttgiation. 


2.  Goidoo,  JJ).,  “AnUyib  and  Hi|ti  Fiequency  Coneetkm  of  the  Ba^Mleaet  Vdodly  Metei  Reapanae,"  DwM 
Taylor  Modd  Baain  Report  2187  (Apr  1966). 

A  conqilete  Batfa^  of  lefemices  it  giren  on  page  SI . 
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APPROACH 


Iidtud  Titoe  and  startup  probiems  wiD  be  soroiounted  by  providing  the  integration 
system  witii  feedback  so  that  steady-state  motion  wiD  be  tq>roduced  at  the  integration 
^stem  on^t  with  die  loss  of  very  low  frequency  cmiqionmits  only.  The  feedback 
tnaagaaent  and  parameters  will  be  diosen  so  tbat  transient  motion  due  to  the  test  wiD 
be  acivratdy  rqirodaced  for  a  finite  time  after  the  b^jnnii^  of  tte  transimit  The 
stea^-state  requirements  necessitate  analysis  of  die  integratkm  system  fitun  the  viewpoint 
of  stability,  and  the  transient  requirements  dictate  die  nature  of  the  feedback  used  to 
stabOize  die  integration  system.  If  die  stabCking  sgnal  being  fed  back  is  die  intqtstor 
output  ddayed  in  time,  die  Integral  of  a  transient  s^al  wQl  be  exact  for  a  time  interval 
measured  from  the  beginiiing  of  the  transient  and  of  duration  equal  to  that  of  die  feed- 
badc  dday.  After  the  duration  of  the  dday,  an  exaggerated  error  is  produced  in  die 
integrated  result  as  die  output  signal  goes  on  to  meet  the  requnements  of  stabiUty. 

Ibe  feedback  dday  may  be  obtamed  in  either  of  two  ways  or  throu^  a  com1»- 
nation  of  both.  The  first  is  to  interpose  between  the  output  and  the  feedback  a  netwmk 
having  capacitors  ndikh  most  be  charged  by  die  output  before  die  output  signal  may  be 
fed  back.  Charging  the  c^acitors  requires  time,  resulting  in  a  dd^  effect  characterized 
by  a  gradual  increase  in  integral  error.  The  second  approach  utilizes  a  tape-recorder 
tran^ort  dday.  As  fidd-test  results  are  usuaQy  recorded  on  magnetic  tape,  an  efficient 
use  of  equipment  would  be  to  use  the  t^  recorder  tran^ort  dday  between  the  record 
and  reproduce  heads  to  provide  die  feedback  dday.  This  second  method  of  obtainiiig 
feedback  dday  introduces  error  only  after  the  duration  of  die  delay.  A  combination  of 
tape  transport  dday  and  capacitor  network  delay  results  in  exact  integration  for  the 
duration  of  the  t^ie  dday,  followed  by  a  gradual  deterioration  in  integral  accuracy. 

The  ideas  devdoped  thus  far  are  qualitative.  Laplace  transform  transfer-function 
tediniques  are  used  to  derive  numerical  values  for  all  parameters  in  the  design  of  the 
stable  intqiration  system.  The  basic  design  technique  is  to  determine  all  design  parameters 
so  that  a  transfer  function  appearing  in  die  mathematical  aiudytis  will  have  die  overall 
diaracteristics  of  a  unit  dday.  This  transfer  fuiKtion  is  called  the  unit-dday  approximation 
and  contains  the  transfer  functions  of  both  the  network  dday  and  the  uqie  dday  as  wdl 
as  the  integrator  feedback  parameters. 
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'ntroug}!  an  invest^tion  of  a  physical  system  described  by  a  staMe  tiansfer 
function  possessii^  a  zero  at  die  c»i|^  it  is  diown  diat  staUe  nndtiple  integration  may 
be  obtained,  havii^  any  desired  d^ree  of  accuracy  for  a  finite  time.  Several  integration 
system  de'ogns  differing  in  die  nature  of  dieir  error  will  be  presented.  Two  system  des^ns 
exhibiting  unidiiectional  error  are  fiilly  treated;  iriieieas,  a  third  system  dedgn  exhibitii^ 
trqilde  error  is  discussed,  and  an  ex«n^  is  ^<il  The  systems  are  investigated,  using  an 
analog  computer,  and  die  practicability  of  die  int^ration  mediod  is  then  demonstrated 
by  urii%  it  to  obtadn  djsfdacement  from  an  accderometer.  The  doutde  integration 
invdhred  is  undostood  to  take  {dace  simultaneously  widi  the  generation  of  the  acceleration. 

THE  STABILIZING  TRANSFER  FUNCTION 

PROBLEM  OF  INSTABILITY 

The  transfo*  function  of  ui  initiaQy  quiescent  linear  system  is  defined  as  the 
ratio  of  the  Laplace  transform  of  die  output  to  the  Laplace  transform  of  die  input.  The 
system  is  stable  if  the  output  due  to  a  finite^tep  input  remains  finite  for  aD  time  and 
is  unstaUe  if  die  output  increases  vrithout  limit 

The  transfer  function  of  the  system  accompliriiing  multiple  integration  is 


vriiere  E2(s)  is  die  transformed  signal  to  be  integrated,  and  E|(s)  is  die  transformed 
result  of  i  integrations.  The  system  described  by  Equation  (1)  is  unstable.  This  instability 
is  die  bask  problem  with  integration. 


Let  the  transfer  function 


=  F(s) 


possess  a  zero  of  order  i  or  greater  at  the  origin  of  the  $-plane.  With  suitable  sdection 
of  F(s),  the  proldem  of  instability  may  be  riiminated  without  si,piificandy  affecting  die 
result  of  the  int^ration  of  a  tranaent  by  constructing  an  integration  system  whkh 
possesses  an  overall  transfer  function.  The  transfer  function  is  the  product  of  correspond* 
ing  sides  of  Equations  (1)  and  (2).  Thus  E|(s),  die  transformed  result  of  exact  int^ration. 
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is  ciiEniittted.  Ihis  dimiiMtioii  leMlts  in 

E;(8)  1 

Ej(t)  s*  (3) 

tdiere  E^(s)  is  t2i6  ^nuisfonned  signal  to  be  integnted,  and  E^(s)  is  die  transfonned  oQ^t 
Hie  system  described  by  F^juation  (3)  is  stabk  became  die  zero  at  tbe  or%in  inbesest  in 
F(s)  cancds  the  pole  at  die  crigin  contributed  by  4-  •  Hietefore  F(s)  is  called  the 
stabilizing  transfor  function. 

Eqoadon  (3)  is  not  to  be  interproted  as  describii^  two  di5*mct  qrstems  in  cascade. 
Equation  (3)  is  die  uransfer  A-’oction  of  one  system  haring  diaiacteristks  of 

two  systems  described  by  Eiputions  (1)  and  (2)  as  ^  diey  were  in  cascade.  Ihe  system 
described  by  Equation  (3)  may  be  regarded  as  per/ornm^  the  intention  descidied 
Equation  (1)  and  then  psssii^  die  result  durougb  the  syct^m  described  by  Equation  (2). 
However,  these  two  operations  are  performed  simoltaneoody,  not  individnsily.  The  form 
sod  parameters  of  Ffs)  win  be  determined  t;*  meet  the  r^irement  that  transients  be 
int^rated  »  accurately  as  possible. 

FORM  AND  PROPERTIES  OF  THE  ST.ABIUZING  TRANSFER  FUNCTION 

Three  dosed-loop  systons,  each  consisting  of  one  or  more  intqpators  with  feed¬ 
back;  an  adder,  a  tape-transport  dday;  and  a  network  ddr^  havn^  poles  only  are  dhown 
in  Figure  1.  The  number  of  integrators  witfi  feed!»ck  shown  for  eadb  system  is  j,  and 
the  differential  order  of  the  system  is  k.  The  transfer  fiiiKtion  of  tiie  tape-tranqiort 
delay  is  exp(-ds),  and  dte  transfer  function  of  the  network  dday  is 

G(s)  = - i - 

k 

I 

a=j 


For  each  syston  shown  in  Figure  1,  the  transfer  function  between  the  irqmt  and 
output  of  tile  adder  is  given  by 


k 


EqCs) 

Ei(s) 


F(s)  = 


Z  V" 


i-i  k 

exp(-«s)  I  J,s-  1  J,s« 

fc=o  *f=j 


<4) 


rriiere  J|t=  1, 

0  is  the  dday  time  of  the  tape  transport,  and 
J,  are  coefficients  to  be  determined. 
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Equation  (4)  is  the  stabilizing  transfer  function.  The  numerator  of  Equation  (4)  has 
a  ZCTO  of  order  j  at  the  ori^  of  die  a^itane  which  cancds  die  ^  in  die  denominator  ofi 
Equatioti  (3)  for  values  of  i  through  j. 


Figure  1  •  The  Analog  Simulation  of  the  Stabilizing  Transfer  Function 


. . . . . . * . * . . . . . .  ii , . I,,.;,  . . . . . . . . 
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Assamii^  fliat  the  system  described  by  ‘Equatioi?  (4)  is  staUe,  tpirikation  of  die 
final  value  theorem  of  the  Laplace  tiansfonn  to  the  transformed  step  reqionse  of  die. 
system  described  by  Equation  (4)  shows  the  characteristics  of:  the  step  re^onse  which  are 
independent  of  the  values  of  die  coefficients  in  Equation  (4).  .^plication  of  diis  <dieorem 
riiows  that  die  step  response  and  its  integrals  through  the  integral  approach  zero  as 

time  grows  infinite.  In  order  for  the  integral  conditions  to  hold,  the  stqi  re^onse  must 
have  appropriate  amounts  of  positive  and  negative  area.  These  areas  are  produced  by  an 
oscillatii^  response.  The  minimum  number  of  oscillations  possiUe  in  die  step  response 
is  that  necessary  to  produce  the  minimum  number  of  positive  and  ne^tive  areas  so  that 
the  integral  conditions  just  described  will  hold.  The  minimum  number  of  (»ciilations  is 
thraefore  determined  by  the  order  of  the  zero  at  the  ori^  of  the  deanrilting  transfer 
function. 

When  the  number  of  oscillations  pre^nt  in  the  step  response  is  equal  to  the 
minimum  determined  by  the  value  of  j,  the  syston  will  be  called  very  staUe.  Vtiien 
the  number  of  oscillations  is  greater  dian  the  minimum  but  finite  in  amplitude,  die  system 
win  be  caded  oscillatory.  Since  unavoidaUe' oscillations  occur  in  die  latter  time  response 
to  an  input  st^,  the  step  may  not  be  reproduced  at  these  latter  times,  and  ati'-nition  . 
must  be  focused  on  stqi  reproduction  in  the  interval  of  time  immediately  after  the  occur¬ 
rence  of  die  step. 

COEFFICIENT  DETERMINATION 

Equation  (4)  may  be  given  in  the  form 


-^=  1  -  Y  II.  s» _ _ 

^iW  n=o-»o  H  k 

^  +  expt  ds»  X 


. . .  . . . 


From  Equation  (5),  it  is  seen  that  the  output  of  the  system  described  by  Equation  (4) 
is  equal  to  its  input  minus  the  responses  of  the  input  and  its  first  j-1  conseaitive  deriv- 
atives  (each  wei||ited  by  )  to  the  system  described  by  tiie  transfer  functiocr 


^(s)  -^o 

£.(8)  i*  k 

+  exp(  0s)£ 


(6) 


If  the  re^nse  of  the  system  described  by  Eqt-ation  (6)  to  the  input  of  the  system 
described  by  Equation  (5)  and  itsj-1  consecutive  derivatives  is  delayed  :n  time,  the  out¬ 
put  of  the  system  described  by  Equation  (5)  will  be  equal  to  its  input  for  the  duration 
of  the  delay.  Therefore  the  J  coefficients  will  be  determined  by  a  criterion  of  approxi¬ 
mation,  which  ^es  the  system  described  by  Equation  (6)  the  essential  character  of  z 
unit  delay.  A  variable  delay  may  be  obtained  by  time  scaling. 

Equation  (6)  wiU  be  called  the  unit-delay  approximation.  It  contains  the  tape 
and  network-d^ay  transfer  furctions  as  we.*l  as  the  integrator  feedbac  k  parameters.  When 
all  quantities  in  Equation  (6)  are  specified,  the  stabilizing  transfer  function  is  completely 
determined.  The  criterion  for  unit-delay  approximation  wfll  now  be  established. 

After  some  manipulations  with  the  MacLaurin  series  of  the  natural  log  of 
Equctim  (6),  Equation  (6)  becomes 


where  the  C  coefneients  are  functions  of  the  )  coefficients  and  6  .  By  makii^  C]  unity 
and  die  next  k-1  consecutive  odd-subscripted  coefficients  zero,  k  J  coefficients  are 
determined,  and  Equation  (7)  may  be  appro.<unated  by 


Ejls) 

Eyis) 


<8) 


Snee  expf-s)  is  the  transform  of  a  unit  delay,  the  unit  impulse  response  of  the  system 
approximated  by  Equation  <81  is  the  inveise  transform  of  the  function  of  s  in  the 
tnackets,  shifted  to  the  ri^t  one  unit  of  time. 
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It  can  be  shown  that  the  double  sided  unrene  tnmsftmn  of 


is  DgCxp 


/'y 


*1 


(9) 


«4t«e  Dgis  diosen  so  that  die  ioea  under  fhe  time  function  i«  unity,  and  die  lemaudog 
D  coeffidents  ate  ctetetmined  itmn  die  C  coefBdents.  Expression  (9)  is  rimilar  to  the 
standard  nomud  curve.  Therefore  die  amt-impolae  re^ioiise  of  the  system  q^Mroximated 
fay  Equation  (8)  mqr  be  regarded  as  a  normal  type  of  curve,  cmitered  about  unit  time. 

As  die  cndo*  of  die  nnit-dday  approximation  inciemec,  the  nomud  type  of  curve  becomes 
higher  aid  narrower,  approadiii^  a  nrnt  impulse  ddayed  one  unit  of  time.  Thus  it  is 
seen  that  use  of  the  criterion  sdected  to  determine  the  J  coefiRdents  gives  die  system 
(Ascribed  by  Equation  (6)  unitdeiay  properties. 

The  J  codiddents  are  readily  e/duated  in  terms  of  n,d,  and  k  if  a  set  of  linear 
equatioits  is  sdved  lAidi  guarantees  dut  conditions  put  on  the  coefficients  of  the  odd 
powers  of  s  in  Equation  <7)  are  met  The  deiivadon  of  this  set  of  linear  equatioiis  foDows. 
Let  the  MacLaarin  expansion  of  die  reci|Wocai  of  Equation  (6)  be  given  by 

j-i  J  k  J  « 

^ --ls«  +  exp(ds)2;  _1s*=2B,s- 

»*«  »=o  (10) 

CbmlHiiii^  Equations  (!0)  and  (7)  yidds 


where  1.  It  may  be  diown  diat  vdien  C]  =  I  md  die  next  k*l  odd^bscripted 
C  coefikiaits  are  zero, 

1  =B, 

f  =  B,-B3 

B3-3B4+3B5 


y  B.^  =  e*P  {  Z  C 
=  1 


e 


I  2al 

- - =y  ...B, 

^o,n  ■=■ 

■  =I,  2,...,k  (12) 

wher?  die  A  coidTkients  are  seen  to  be  the  Sessd  Po)yno'.iid  co^doits-'  and  are  ^tn  by 

,  _  2»(2^)! 

Aofl - - - 

2Pa!(^)! 

Refcicace  3  -  KnI,  HX.  a»d  Onte  Fiiak,  "A  N«w  Oik  of  OtlkotOMl  PoIjnMWiWK  Hw  Scael  PoiyfiUs." 
llainctiOM  of  Ike  AaeiicM  Milhcmalica  Society,  VoL  SS.  No.  I.  pp.  100-107  (Jw  1949). 

A  cowphte  liM  or  lefcfMKU  b  fhea  om  ptf/e  SI. 
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(13) 


'Ihe  B  coeffideats  say  be  jgKVi  n  temis  ctf  the  J  codfldmti  »d  6  by 


B,*  -^’  o<u<j-l  I 

\ 

y - :  '7“»»>v  \ 

::,  <•■"' ; 


SabfititiiticMi  <rf  EqwtioBs  (13)  in  EquatkNS  (12)  gives 
1  2a-l 


— -I  (*!)•*■  ..I- 

A»,*  •*m  *( 


.j<S 


•  *  i . H 

k  2«-l 


-J-=S  I  M)-"  Aj.., * 

A,ji  .11  ».« 

3  }>▼ 


»•«  »»i  »*• 

•*i*  -  •  -f  ^ 


(14) 


The  ripyinn  of  (14)  gives  k  imear  ansaltaiieoai  eqaalkHiB  in  k  ndinown  J 

coefficient  ntios.  After  ^ledfytng  a  ptftkaitf  vatee  of  6,  lhae  eq^tions  may  be  sdvei 
for  tfie  J  co^Bcknts  cone^onding  to  ^t  valne  of  g. 


MAC  LAURIN  SYSTEM 


A  systea  des^ned  frewa  a  nnit-dday  igqHoxijaation  hai^  zero  9  wffl  be  cdled  a 
MacLooiin  syatem.  Sokitkm  of  Equations  (14)  for  zero  ^cs  t*ie  J  coefficients  forMacLmuin 
systems  dneedy.  They  m 


t 

■*« 


■  A,.k 


2*(2k-m)! 

2kii!(k-.)! 


(IS) 


Table  1  lists  die  J  coefficients  few  MacLaorin  systems  widi  k  ruiging  froin  1  to  10. 


TAPE-DELAY  SYSTEM 

A  systmn  designed  from  a  onit-delay  approxiniation  havh^  9  greater  due  zero 
artt  be  c^ed  a  tafie-dday  system.  The  sohition  of  Eonations  (!4)  widt  9  varied  diows 
dut  if  (k^)  is  odd.  Jg/Jo  P*®®  dwoof^  zero,  gmm;  native  as  g  is  tscieued.  The 
system  is  very  stdile  nnti  iiBtabffity  occcr  'Sith  die  cha^  m  of  ig/io-  Therefore 
$  nuy  not  be  fnrdieT  increased.  For  the  ’'slue  of  k  md  9  iwoAicing  zero  Jg/Jc-  die 
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remaiiiuig  J  co^cient  ratios  are  to  the  J  coefSkient  ratios  cakuiated  fcH'  the  next 
lower  order  and  the  sune  valo3  of  6.  Thus  two  adjacent  orders  ^e  die  same  unit  delay 
ipproxhnation  for  a  certain  v^e  of  d .  As  0  is  increased  fr<mi  diis  valce^  die  Iowct  of 
the  two  orders  becomes  osdDatory.  Therefore  die  value  of  0  prodocine  zero 
(k-j)  odd  k  die  opthmuD  value  of  6.  Ffr  9  opthnnm  Table  2  Usts  6  and  die  corre^nd* 
ing  J  coeCBcients  for  the  first  duee  values  of  j  with  k  rangii^  firmn  j  to  10. 

STEP  RE9*ONSES 


Now  diat  the  J  coefficients  have  been  deteimined,  and  die  stalnUzIrig 
function  of  Equations  (4)  and  (5)  is  cimipietdy  defined,  die  diaractimstks  of  the  stabi¬ 
lizing  transfer  function  wQl  be  studied  throng,  analog  shnulation. 

Figure  1  shows  die  analog  Emulation  of  die  stabilizing  transfer  faartion  for 
the  first  three  values  of  j.  By  recording  d’.e  output  of  the  adder,  or  magnetic  tape,  and 
by  ^ultaneously  placing  the  output  back  into  the  input  of  the  simolation  of  G<s),  a 
^nal  ddaqr  equal  to  the  time  required  for  the  tape  to  pass  from  die  record  head  to  the 
reproduce  head  wiD  be  obtained.  Fordier  dday  in  die  MacLaurin  sense  is  obtaiiied  with 
the  shnulation  of  G(s).  For  convenieiKe  in  simulation,  G<s)  be  factored  anJ  given 
ni  die  form 


G(s)= 


U3 


2 


(16) 


where  die  brackets  indicate  fractions  are  dropped.  Table  3  gives  the  parameters  of 
Equation  (16)  for  the  MacLaurin  systons  of  Table  1  for  the  first  three  values  of  j. 
Table  4  gives  die  parameters  of  Equation  ( 16)  for  the  tapenlday  systems  of  Table  2. 
Rgures  2  through  4  are  step  responses  of  die  MacLaurin  system  obtained  by  analog 
simulation  for  j  =  1,2,  and  3,  reqiectivdy,  with  k  rai^h^  from  j  to  10.  F^ures  5 
through  7  are  the  initial  portions  of  Figures  2  throu^  4,  respectively,  widi  scales 
expanded.  Figures  8  through  10  are  stq[>  responses,  obtained  by  analog  simulation  of 
the  tape-dday  systems  given  in  Table  2  for  j=I,  2,  and  3,  reflectively.  Figures  1 1 
dirou^  13  are  the  initial  portions  of  Rgures  8  throng  10  with  scdes  expanded. 
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3*00000 


5*23548 


7*47150 


9*71062 


11*95128 


0*77459 

1 

3*87298 

0*53195 

4*47833 

0*90889 

6*20638 

0*37533 

4*95309 

0*75293 

6*99320 

0*28101 

5*33374 

0*95157 

8*49114 

0*22662 

5*62469 

0*60353 

7*69695 

0*85382 

9*36765 

0*19545 

5*82274 

0*47852 

8*37620 

0*97008 

10*75576 

0*17485 

5*94309 

0*74550 

10*15438 

0*38222 

9*04371 

0*90399 

11*68398 

0*15848 

6*01739 

0*64199 

10*88054 

0*31245 

9*67548 

0*97971 

13*01182 

0*14466 

6*06790 

0*82407 

12*52995 

0*54853 

11*56306 

0*26397 

10*25041 
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TABLE  3  •  (Continued) 


j 

k 

L 

n 

Tn 

2 

3 

6.00000 

2 

4 

00 

1 

0.74939 

2 

"  5 

8.45665 

1 

0.4o424 

2 

6 

oo 

1 

0.90217 

4* 

2 

0.27419 

2 

7 

10.85906 

1 

0.73660 

■* 

2 

0.14983 

'2 

8 

eo 

0.94892 

■ 

0 

0.06739 

0.57991 

2 

9 

13*23045 

0.84638 

/ 

• 

2 

0.01137 

3 

0.44868 

2 

10 

00 

1 

0.96873 

2 

-0*02807 

3 

0.73319 

4 

0*34558 

2 

11 

15*98141 

1 

0*89990 

2 

-0*05716 

3 

0*62528 

4 

0*26137 

2 

12 

OO 

1 

0.97892 

2 

-0.07961 

3 

0*81690 

4 

0*19806 

5 

1  0*52815 

- - 1 


ft97dS29  " 
7.0%7jt5 

9*«2i44 
7*2fS97 
lOa  11.039 
7«49St9 
ll«t9864 
7«69384 
10 *87960 
12*73167 
7.77160 
ll»lt927 
:*«S021? 
7*65238 
13 *64594 
11*66218 
15*21222 
7*90414 
14*08678 
12*10229 
16*66992 
7*93525 
16*01552 
12*50480 
14*67993 
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TABIf  3  •  (Cimtinaed) 


k 

a 

ta 

10*00000 

oo 

1 

0*73192 

10*24699<' 

12*25m 

1 

0*43124 

10*14029 

40 

1 

0*89844 

13*07732 

2 

0*22339 

10*06909 

14*610iS 

1 

0*72494 

13*99983 

2 

0*08347 

10*03011 

04 

1 

0*94738 

19*99397 

2 

-0*01240 

10*01410 

3 

0*94398 

13*92430 

10.9S404 

0 

1 

0*84183 

14*31838 

2 

-0*08048 

10*00838 

3 

0*42442 

14*24218 

00 

1 

0*94791 

18*03149 

2 

-0*12984 

9*99882 

3 

0*72927 

14*90923 

4 

0*31914 

14*93723 

19* 39499 

1 

0989732 

18*87941 

2 

-0*14684 

9*98431 

3 

0*61404 

17*42477 

4 

0*22999 

14*81909 

OO 

1 

0*97842 

20*43492 

2 

-0*19941 

9*96481 

3 

0*31221 

19*40343 

4 

0*19497 

19*07403 

- 

LL 

9*913a2 

17*89909 

TABLE  4  •  NETWORK  DELAY  PARAMETERS  FOR  TAPE4>QLAY  SYSTEMS 

b*  •• 


0«370a0 

0*3»000 


OtiUio 


0«33110 


0«2«470 

0«2S470 


0*26100 


0*26630 


0* 17640 
C* 19640 


0*2096C 


0*26123 

0*14074 

0*46643 

0*09962 

0*24364 

0*66021 

0*07662 

0*16293 

0*36320 

0*66166 

0*16144 

•0*04430 

0*41637 

•0*12069 

0*17171 

0*66014 

•0*16763 

0*06266 

0*32091 

0*64662 

0*10727 

•0*13246 

0*39693 

•0*23066 

0*13032 

0*64732 


6*64169 

6*06491 

11*07396 

6*16060 

11*33769 

16*46313 

6*22176 

11*66626 

16*09377 

22*12326 

6*32236 

6*01464 

14*46764 

7*96906 

13*60260 

20*61346 

7*96906 

13*94196 

19*21329 

26*66772 

11*41247 

10*32666 

16*46666 

10*01661 

16*72666 

26*20036 
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-  Step  Reeponsce  of  MacLaurin  Syttcma  witir  j-2 


Figure 


24 


26 


dais  - 
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Figure  13  •  Figure  tO  with  Scales  Expanded 


A  conqiaiison  of  tiie  step  responses  ^ows  the  manner  in  which  oscillations 
between  Ae  initiall>  flat  portion  of  step  repiesnstation  and  the  eventual  zero  value  increase 
with  increasing  order  of  the  system  transfer-function  zero  at  the  origin.  It  has  previoudy 
been  sbown  that  these  oscillations  must  be  present  so  fliat  the  fhst  (j-llcorsecutive  integrals 
of  the  stq)  response  will  approach  zero  with  increasing  thne.  The  figures  show  that  by 
increasit^  tfie  order  of  the  sy'atem,  tiw  duration  for  udtich  die  output  follows  die  stq> 
input  is  increased.  It  is  apparent  from  the  figures  that  for  a  given  value  of  j,  improvement 
occurs  at  a  decreaai^  rate  as  order  is  increased.  Bov>ever.  for  a  given  deviation  fnnn  the 
euct  step,  percentage  inqrrovemerrt  betwemi  two  given  consecutive  values  of  k  increases 
with  the  vahte  of  j. 

THE  STAfilUZING-SYSTEM  RESONANCE 

From  a  study  of  die  unit-step  re^nses-.  it  is  seen  that  die  amplitudes  of  die 
excursions 'die  signal  takes  grow  with  order  k  and  order  of  system  zero  j,  after  die 
initial  intervals  of  step  transmission.  A  transient  response  of  this  nature  implies  a  resonant 
bnildop  of  oat{N!t  at  certain  input  frequencies.  The  stalnliziiig-systeni  output  may  be  a 
steady-state  signal  such  as  that  produced  by  srater  srave  motion  of  a  diip  at  sea.  To  avoid 
resonant  distortion  of  the  output,  die  stabilizing  system  must  be  des^ned  to  put  die 
resonance  region  at  frequencies  lower  than  dtose  of  interest  in  die  motion.  To  provide 
information  about  resonant  bufldnp,  die  steady-state  sinusoidal  fiequoicy  responses  of  die 
MacLaurin  stalnlizii^  syston  for  j=l,  2,  and  3  are  given,  respectivdy,  in  Flgnres  14 
dirough  16  for  crdeis  j  'through  10- 

TIME  SCALING 

The  duration  for  nliich  die  output  fdUows  the  step  input  may  be  varied  by  tune 
scaling  the  stabOizii^  system.  Time  scaling  is  accomplished  by  substitutiiq;  Tp  for  r  in 
the  describing  transfer  function  and  amulating  die  result  The  dday  tune  of  die  resulting 
system  is  the  time-scale  factor  T.  By  molti{riyiiig  die  time  axis  of  die  step  responses  by 
T  and  by  dhriding  die  angular  frequency  axis  of  the  fieqoaicy  responses  by  T,  the 
rmponses  of  tone-scaled  ^sterns  are  obtainedt 
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Fiffire  k4  -  Frequency  Responses  of  MacLaurin  Systems  with  j=l 
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Figure  15  •  Freuusncy  RaqMnsas  of  MacLaurin  Systaim  wHh  j«2 
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Fi^re  16  -  Frequency  Responses  of  l^Lanrin  Sysfems  with  j^3 


THE  STABLE  INTEGRATION  SYSTEM 


ANALOG  SIMLiLATlON 


With  icference  to  Equations  (3)  and  (4),  it  is  seen  that  the  transfer  function  of 
die  tinte'^c^ed  stable  alteration  system  is 


£o(p)  _  t 
t^(P)  P* 


I 

_ _ _ 

M  k 

exp(-«Tp)  5;  J.(Tp)«+X  J.(Tp)« 


where  Tp  is  substituted  for  s  in  die  stabilizing  transfer  function.  Figure  17  shows  the 
aiming  ^ulation  of  Equation  (17)  for  die  first  three  values  of  i  and  j.  From  Rgure  17 
it  is  seen  diat  the  value  of  i  determines  the  interator  into  whidi  die  input  is  fed,  i  being 
less  d>an  or  equal  to  j.  The  value  of  (j*i)  determines  the  final  vdue  of  die  output  for 
a  given  constant,  i  imp,  parabolic,  etc,  input 

When  (j*i)  =  0.  die  final  value  of  the  output  for  a  constant  input  is  a  constant. 
When  0'*i)  ~  L  the  final  value  of  die  output  for  a  constant  input  is  zero.  If  the  transient 
to  be  multqdy  integrated  has  an  additive  constant,  die  output  of  the  int^pration  systnn 
having  (j-i)  ~  1  wiQ  be  ind^iendent  of  this  constant. 

APPLICATION  TO  MOTION  MEASUREMENT 


The  pracficability  of  die  staUe  muld{de  integration  system  tray  be  demonstrated 
by  uang  the  system  to  obtain  dsidscement  by  double  integration  of  an  acederometer 
output.  The  accuracy  of  die  int^ation  system  may  be  verified  by  comparing  the  result 
to  a  direct  measurement  of  die  displacement.  Figure  18  is  a  picture  of  an  acederometer 
attached  to  a  linear  potratiometer  and  constrained  to  ncove  vridiout  rotation  alor^  a 
piide.  Ihe  dhqplacentent  of  the  acederometer  as  it  is  moved  along  the  gohle  is  meKorsd 
diiecdy  by  die  potentiometer  and  inditeedy  by  double  intqpation  of  the  ecederometer 
output  usii^  an  integration  system.  The  integration  systems  demo.8trated  will  have 
i=2  and  j«3.  Since  (j*i)  =  1,  the  output  is  independent  of  acederometer  zao  shifts  and 
attitude  char^  taldng  {dace  after  the  dectronics  start  but  wdl  before  the  initiation  of 
the  motion. 
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Figure  18  -  Apparatus  for  tl>e  Direct  and  Indirect  Measurement  of  Displacement 


MAC  LAURIN  SYSTEM 


Figure  19  is  tiie  analog  simulation  diagnm  of  the  ntntti  order  MacLaniin  qrstem, 
time  scakd  for  T^27.4.  Figure  20  shows  graphs  of  disphonnent  obtained  simoltaneoudy 
from  both  tile  potentiometer  and  the  doubly  intqrated  accderomet^  as  die  accderometer 
is  moved  along  the  guide.  A  cmnparison  of  tiiese  two  gn^iUi  diows  the  manner  in  whkdi 
error  accumulates  in  tiie  douUy  integrated  accdennneter  output  whoi  a  MacLaurin  qrstem 
is  used.  Theoieticaily,  for  a  AfacLanrin  system,  error  accumulates  starth^  at  the  b^nnii^ 
of  die  motimi.  However  tfie  error  accumulates  very  graduMly  and  only  becmnes  signifi¬ 
cant  at  later  times.  For  a  step  diqdacement,  tiie  system  being  demonstrated  accunmlates 
an  emw  of  ^roximately  0.75  percent  at  S  sec. 

TAPE-DELAY  SYSTEM 

Fq^re  21  is  the  analog  simulation  diagram  of  tiie  ninth  order  tape-dday  system 
with  0^.2095  sec  and  T»27.4.  Figure  22  diows  graphs  of  diqdacemoit  obtained  amnl- 
taneously  from  both  tiie  potentimneter  and  the  accelerometer  double  integraL  From  these 
graidis  it  is  seen  tiut  dnoarmble  error  does  not  accumulate  until  about  10  sec  after  initi¬ 
ation  of  motion.  This  10  sec  of  double  integration  is  the  sum  of  5.74  sec  of  exact  double 
integration,  controlled  by  the  tape-transport  dday,  and  4.26  sec  of  double  integration 
with  gradually  increasing  error,  contndled  by  tiie  network  dday. 
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Figura  19  •  Tht  MacLaurin  Doubia  Intagration  Sysiam 
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Figure  20  -  Direct  and  Indirect  Displacement  Measurements,  MacLaurin  System 


Figura  22  •  Direct  and  indiract  Diiplacamant  Maasuramant,  Tapa-Dalay  Syctam 


ALTERNATIVE  SYSTEMS 


FADE'  EXPANSION  OF  exp(-ds) 


In  the  preceding  sj^tems,  die  physical  device  having  die  transfer  function  exp!-ds) 
was  die  t^  tran^ort  A  completely  dectronic  approximadon  of  exp(-ds)  may  be  obtained 
by  using  die  Fade'  ex{mnsion^  of  exp(-ds),  having  numerator  and  denominator  polymnnids 
of  equal  degree.  This  Fade'  expansion  is  given  by 


exp(-ds)  « 


(18) 


where  q  is  die  degree  of  the  numerator  and  denominator  polynomials  of  the  expansion, 
and  the  A^  ^  coefficients  are  the  same  as  appear  Equations  (12).  .Vhen  the  rq^t-hand 
side  of  Equation  (18)  is  given  as  e  raised  to  a  power  series  in  s,  q-l  consecutive  odd 
powers  of  s  after  the  first  are  zero  as  in  Equation  (8).  Therefore  the  Fade'  expansion  of 
exp(-ds),  having  numerator  and  denominator  polynomials  of  equal  degree  satisfies  die 
criterion  of  dday  approximation  previously  establidied. 


RIFFLE  SYSTEMS 


When  design  data  for  tape-dday  systems  are  used  with  the  simulation  of 
Equation  (18)  teplacing  the  ti^ie-transport  dday,  die  resulting  integration  system  has  an 
output  widi  oscillating  error  in  the  initial  time  interval.  By  increasing  q,  the  amfditude 
of  die  error  <»cQlations  are  reduced  and  may  easily  be  made  ne^igiUe.  Since  die  output 
recording  device  is  not  used  as  part  of  the  integration  system  when  the  simulation  of  die 
Fade'  expansion  replaces  die  tiqie  transport,  an  additional  feedback  may  be  made  at  die 
input  of  the  exp(-Os)  simulation,  resulting  in  an  increase  in  the  time  interval  of  integra¬ 
tion.  A  syston,  different  from  the  tape-dday  systan  already  treated,  results  from  adding 
die  extra  feedback;  thus,  a  new  unit-delay  approximation  must  be  derived,  and  its  nara- 
meters  must  be  determined  using  Equations  (12). 

Rohranca  ••  -  rrunl.  John  G.,  “AutoiMtic  rwiawLfc  Contra)  SyfMm  Symtiwii,"  McGraw-HiN  Boole  Co,  Ine, 

Him  Yoifc,  ItY,  (1965)  p.548. 
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As  an  e:;ample  of  an  integration  system  having  ripple  ernH-,  the  syston  havh^ 
the  additional  feedback  will  be  treated  uidng  an  output  agnal  dday  consisting  of  the 
sinnilation  of  the  fourth  degree  Fade'  expansion  of  exp(-ds),  foUowed  by  a  second 
order  networic  delay  having  pdes  only.  Figure  23  is  die  analog  ^uladon  dh^iram  of 
die  douUe  integration  system  with  j=3.  The  unit*dday  i^iproximation  analogous  to 
Equation  (6)  for  this  system  is 

Ms)  ^  _ _ Jo _ 

£4(5)  Jo+J,s+J2s2+s3|j3  +^i£lJ4i!^j.exp(fls)j 

vdiere 

Jo  =  8.247 
J|  =  8.247 
J2  =  3.915 
^3  =  .16649 
J4  =  129.14 
J5  *  —.8959  and 
9  =  .3497 


are  determined  from  Equations  (12).  When  a  dday  of  6  seconds  is  approximated  by 
simulation  of 


exp(-ds) 


105- 105 1 

[■ff  «( 

-10 

105  +  105 

(t»)'  « 

die  output  step  in  the  initial  interval  has  an  oscillating  error  of  only  0.25  percent  as 
shown  in  Figure  24.  If  a  hi^er  degree  Fade'  expansion  for  the  0  dday  had  been  simu¬ 
lated,  die  number  of  oscillations  in  the  interval  of  step  representation  would  have  increased 
and  dieir  amplih'de  would  have  decrsued. 

Figure  25  shows  the  results  of  an  accderometer  double  int^ration  test  using  the 
ripide  system  example  shown  in  Figure  23  with  T=27.4.  From  Figure  25  it  b  seen  diat 
ripple  error  is  ne^igiMe  for  about  9  sec  aftm-  die  initiation  of  the  motion.  After  die 
lapse  of  9  sec,  the  response  goes  on  to  meet  the  requirements  of  stabflity  in  the  same 
manner  as  in  the  previous  systems  discussed. 
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Displacenwnt  Meaniremento,  Ripple  System 


. . . . . . 


DISCUSSION 


The  success  adiieved  ai  obtsiniiig  diqiliceinent  by  doubly  integntay  an  accder-- 
ometer  output  diows  that  the  staUe  int^ntioii  method  works,  nsii^  die  feedback  ainuige> 
ments  tried.  Using  the  taperecorder  transport  delqr  in  tire  int^ratkm  system  results  in 
a  dispiacsnenl  record  having  no  error  fo.'  die  duration  of  the  traruport  dday;  however, 
usii^  die  network  dday  alone  ^es  a  satisfactmy  record,  even  diou|jh  error  budds  up 
from  die  begiitriing  of  die  motion.  This  buildnp  of  error  is’  very  gadiud  at  first  and  is 
kept  widkin  acceptable  limits  by  rairing  system  order  and  Iqr  time  scaling. 

When  qnasi-periodic  transtiocer  motion  is  present  die  integration  sy|deiii  resonance 
nmst  be  considered  in  the  system  des^  This  considefatioa  affects  the  selection  of 
system  order  tnd  tnne  scdh^  and,  dierefore,  system  error.  The  resonance  of  the  system 
is  increBaed  when  die  error  is  reduced  by  nush^  system  order,  and  the  settling  time  of 
die  system  is  increased  wbrn  die  error  b  reduced  by  time  scaling.  No  matte  how  error 
comprosiiaes  are  made,  die  requircmentii  of  stability  intrcdnce  long  time  enor. 

li;dusion  of  a  tape^tnuE^liort  dday  in  the  'eedback  arrangement  may  be  regarded 
as  a  madiematicd  design  technique,  ossurit^  zero  error  for  the  duration  die  dday. 

The  actud  tran^iort  dday  does  not  have  to  be  pn^sent  It  can  be  approximated  dectron- 
ically,  resulting  in  an  integration  s}rstem  having  osciHadr^  error. 

Three  of  tiie  nudn  factors  involved  in  sdectii^  a  particular  int^mtkm  system 
design  ne  quantity  of  dectronics  needed  for  fabrication,  rdiabflity,  and  system  error.  Of 
the  diree  analyzed  die  tapedday  system  requires  the  most  equipment  and  the  most 
maintenance.  The  tape  transit  b  inherendy  less  rdiaUe  dian  completely  dectronk 
devices.  However,  die  tape-dday  system  has  die  minimum  oror  of  the  systens  dbcassed 
since  its  prodnces  no  error  for  die  duration  of  the  transport  delay.  The  MacLaurin 
qrstem  requires  die  least  equipment  and  b  die  most  rdiaUe.  The  dtsadvantige  b  diat 
it  has  die  smallest  ratio  of  useful-to^tdii^  time  of  the  systems  dbcuced.  The  ripfde 
systems  dbcnssed  are  completdy  dectronk  versons  ot  sysiems  vdiidi  could  also  employ 
a  tape-transport  delay.  The  amount  of  rip|de  ent>r  deprods  on  the  feedbadc  anai^emeRt 
as  wed  as  die  quantity  of  dectronics  used  to  sonulate  die  tapr-transport  dday.  When 
ssnulation  of  the  tapf-tranqiort  dday  b  fdlowed  by  a  second  order  fSte,  the  ripple 
eiror  b  attenuated  to  an  acceptable  levd.  The  feedback  srnngemmt  of  lipfde  systems 
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must  be  built  to  clos^  tMennces  dun  ue  requiied  for  tq>e-delay  and  MacLanrin  systems. 
However,  die  laiger  ratio  of  oseftil>to-settlbig  time  obtained  widi  iqiide  qratems  may  be 
worth  die  additional  effort  required  to  achieve  those  tolerances. 

SUMMARY 

Using  ddayed  feedback  to  stabilize  electronic  integiatcHS,  aQows  the  io^iators 
to  produce  an  accurate  mult^e  integral  of  a  transiait  eiecvricai  sigoal  widiout  it  beii^ 
necessary  to  start  die  integrators  at  the  begituiiiig  of  the  s^nal.  Altfaou^  transient  dec- 
trkal  signals  of  any  origin  may  be  integrated,  using  the  integration  method  presented,  the 
mediod  is  particolarly  iqqdicaUe  to  the  integration  of  electromeclianicai  tra]is(h*ceroutpot 
sigtuds  simultaneously  with  dieir  generation.  The  t»oad  range  of  integration  system  design 
dat<  ^en  enables  the  designer  of  measuremmit  ^tems  to  matdi  the  most  suitable  mte- 
gration  system  widi  the  transducer  sensing  die  quantity  vdiose  int^ral  is  dedred.  The 
most  apparent  application  of  the  integration  method  lies  in  the  area  of  absolute  motion 
measumnent  using  inertial  instrumaits.  The  des^  procedures  and  techniques  presented 
herein  are  direedy  applicable  to  measurement  systems  designed  to  be  used  in  exidosion 
experiments  and  provide  sdudons  to  critical  proUems  associated  wld;  designing  such 
systems. 
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